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Abstract We study the fine splitting in charmonium spectrum in the quark model with the channel coupling 
effect, including DD, DD* , D* D* and D s D a , D S D*, D*D* channels. The interaction for channel coupling is 
constructed from the current-current Lagrangian related to the color confinement and the one-gluon exchange 
potentials. By adopting the massive gluon propagator from the lattice calculation in the nonperturbative 
region, the coupling interaction is further simplified to the four-fermion interaction. The numerical calculation 
still prefers the assignment 1 ++ of X(3872). 
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1 Introduction 

A series of hidden charm states, the so-called A, 
Y , Z, have been discovered and confirmed by the ex- 
periments since 2003. The nature of these narrow res- 
onances has attracted much attention, because their 
properties are not consistent with the prediction of 
the quark model. 

The typical A (3872) state, which was discovered 
in 2003 by the Belle Collaboration [l[ and subse- 
quently confirmed by the CDF Collaboration Q and 
BABAR Collaboration etc., is now listed with 
M x = 3872.2 ±0.8MeV, T x = 3.0l H ± 0.9MeV in 
PDG0. Its quantum numbers were inferred J PC — 
1 ++ or 2 h . The corresponding charmonium candi- 
date in the quark model is 2 3 P 1 or I 1 Z3 2 respectively. 

The mass of the 2 3 P 1 state in the quark model is 
~ 100 MeV above M x . However, the channel cou- 
pling effects by the creation of open charmed meson 
pairs can produce significant mass shift to the bare 
charmonium spectrum. In Ref. Q , only the fine split- 
ting in the mass shift induced by open-charm states 
is considered. In Refs. 0, 0], the whole mass shift 
is considered to lower the bare mass of the excited 
charmonium state. The mass shift can be also hand- 
ily treated by introducing screened potential into the 
quark model[8j. 



The proximity of the X(3872) to DD* threshold 
implies that the cusp scenario may be important 
The cusp can be calculated from channel coupling 
and the result is in qualitative agreement with ex- 
periment [nj. The observed but Okubo-Zweig-Iizuka 
(OZI) forbidden decay channel pJ/ip is also consid- 



ered in Ref. 



Recently, a study of the n + ir ir° mass distribution 
from the A (3872) decay by the BABAR Collabora- 
tion favors the negative parity assignment 2 h 12j. 



However, the mass of the corresponding charmonium 
state 1 1 £> 2 in the quark model is ~ 100 MeV below 
M x . Since the ip(3770) is assigned to 1 3 D 1 in the 
quark model, the assignment 2 h seems to conflict 
with the small fine splitting in cc ID multiplet from 
the quark model calculation 13 1. 

The mechanism of channel coupling is the same 
as strong decay's. The simplest decay model is the 
so-called 3 P model based on the flux-tube-breaking 
model 1J, [l5| . Another model is the Cornel model 
which tries to relate the pair-creation interaction to 
the potential in the quark model IB, 12 1. The Cornel 
model assumes the Lorentz vector confinement so the 
total vector potential is 



V{r) = -- + - 2 . 

r cr 



(1) 



Thus in the Cornel model the decay amplitude from 
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the one-gluon exchange and that from the confine- 
ment add destructively. A similar calculation but us- 
ing the Lorentz scalar confinement shows that the 
decay amplitude from the scalar linear confinement 
is too large 

The lattice calculation shows that the gluon prop- 
agator is quite different in the nonperturbative re- 
gion. The gluon may get a mass about 600 ~ 1000 
MeV 19l421|. A non- vanishing gluon mass is used 



in the phenomelogical calculation of the diffractive 
scattering 22] and radiative decays of the J/ ip and 

In this work, we will consider the fine splitting in- 
duced by channel coupling with open-charm states, 
including DD, DD% D*D* and D S D S , D S D*, DID*. 
Following the Cornel model, we will construct the 
model pair-creation interaction from the potential in 
the quark model, i.e. the scalar confinement plus the 
vector one-gluon exchange. With the assumption of 
the massive gluon propagator in the pair-creation pro- 
cess, we will obtain a simple effective four-fermion in- 
teraction which is quite similar to the case of weak 
interaction. In Sec. 2, we will introduce the channel 
coupling model. In Sec. 3, the numerical analysis is 
performed. Finally, we will give a brief summary. 

2 The channel coupling model 

In the simplest version of channel coupling model 
[7], the hadronic state is assumed to be represented 

by 

!*->=( ° M V (2) 

where the bare state \ip a ) is coupled to several meson- 
meson channels |Afi(i)Af 2 («)}. The system Hamilto- 
nian reads 

"' ' V (3) 



H 



V H 



Ml AI 2 



where H c is the meson Hamiltonian of the quark 
model, with 

H c \iP a ) = M a \iP <x ). (4) 

In this work, H MlM2 includes only the free meson 
Hamiltonian, so 



The Hamiltonian in the non-relativistic quark po- 
tential model can always be written as [7J 



H c — H Q + H sd , 



(6) 



where H and H 3d are the spin-independent and spin- 
dependent parts respectively. The spin-independent 
part reads 

P 2 



2[i 



+ V{r) + C, 



(7) 



/i is the reduced mass. The potential V(r) is usually 
taken to be a sum of the linear confinement plus the 
one-gluon exchange Coulomb potential: 



V(r) 



4a s 

■ ar 

3 r 



(8) 



H sd includes spin-spin, spin-orbit and tensor force: 

H sd = • S 2 + V LS (r)L -S + V T (r)T, (9) 

which determines the fine splitting in the spectrum. 

The off-diagonal interaction V is responsible for 
channel coupling. It depends on the pair-creation 
mechanism of the specific hadron decay model. The 
3 P () model 14, li[ and the Cornel model 

0,113 

arc 

two popular decay models. 

To describe the creation of a light-quark pair in 
the quark model, a plausible approach is to consider 
the quantum field expression of the quark potential 
V(r). In the Cornell model, the quark pot ential is 
replaced by an instantaneous interaction [id . ET 



Hl = \j ^ x ^y-p^ x )\ v ( x -y)p^y)- 



where 



(10) 



(11) 



Hm x m 2 — H Ml + H M2 . 



(5) 



is the quark color-charge-density operator, and ip( x ) 
is the quark field operator. As the spin splitting in 
charmonium spectrum and the lattice gauge calcu- 
lation indicate that the confinement current should 
be the Lorentz scalar, in Ref. [3] the instantaneous 
interaction is replaced by the scalar confinement in- 
teraction plus the vector one-gluon exchange. 

Following the Cornel model, here we will model 
the pair-creation from the quark model. We first as- 
sume the nonlocal current-current action of the quark 
interaction [241 ]: 
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1 f - 1 - 1 

If - 1 - 1 

~2 / d4 z d4 #( x )2 A aV^)s(^y)v^y)-A ?/%)- 



(12) 



The vector kernel G is obtained from the one- 
gluon propagator. In the momentum space 

Ana, 



G(q 2 



(13) 



The scalar kernel S{x — y) is obtained from the linear 
confinement 

S^) = -"f. (14) 

The lattice calculation shows that the behavior of 
the gluon propagator is quite different in the nonper- 
turbative region. The gluon may get a mass about 
600 - 1000 MeV With the gluon getting a 

mass in the nonperturbative region, we can make the 
non-relativistic approximation q 2 — > q 2 ~m 2 « —m 2 in 
the quark-antiquark pair-creation process. Thus 



D»„(q 2 ): 
D(q 2 ): 



ni- 



3 



(15) 
(16) 



Then the channel coupling interaction is simplified to 
the four-fermion interaction 

1 Aita f - 1 - 1 

2 m 2 J 2 2 

+ \^ J d 3 x4,{x) l -\ a ^)^)\K^{x) (17) 

Once we calculate the transition amplitudes 

f i (p) = (ip a \V\M 1 (i)M 2 (i)), (18) 

where p is the relative momentum between Mi and 
M 2 , the mass shifts are given by 



ff (M)=X>(M), 

i 

fi(p)fi(p) 



i(M) = 



2M i 



(19) 



-d 3 p, (20) 



where m-n and m i2 are the masses of M 1 (i) and M 2 (i) 
mesons, ^ is their reduced mass. 

To calculate the coupling matrix element, we will 
use the simple harmonics oscillator (SHO) wave func- 
tions as usual. The partial-wave amplitude f la can be 
expressed as 



F ls (p), 
(21) 



where j3 B = (3c, m c is the mass of charm quark, m q 
is the mass of light quarks (u, d, or s). F ls (p) is a 
polynomial of p which depends on the specific channel 
(the formulas are collected in Appendix). 

Our calculation is basically non-relativistic. But 
the exponential factor in the obtained partial-wave 
amplitude Eq. (f2~Tj) is obviously not enough to cut of 
the high momentum contribution. We will make an 
additional cutoff to the momentum integration. The 
mass shift is then replaced by 



9i(M): 



Up) Up) 



(run 


+ m i2 


+ PL) 


— M 






^ 2« J 





exp(-p 2 /A 2 )d 3 p, 
(22) 

where A is the cutoff parameter. 

Since the channel coupling calculation is essen- 
tially the virtual charmed meson loop calculation, the 
quark potential in the quark model should be renor- 
malized 0] ■ The renormalization process can be out- 
lined as follows. The full Hamiltonian is divided into 



Hfuu =H C + AH. 



(23) 



H c is the original quark model Hamiltonian. Its spec- 
trum is given by 

M nslj =M nl +(V HF ) (S 1 -S 2 )+{V LS ) (L-S) + (V T ) (T), 

(24) 

where M nl is the centroid of nl multiplet which is 
obtained from the spin-independent Hamiltonian H 
and the remaining terms give the fine splitting. (T) 
is the expectation value of the tensor operator, 



i i+i 

6 2i-l 



(T) = 



(25) 



■ - — — 7 = I + 1 

6 2i+3 J t + i ' 



where the total spin s = 1. AH is the cancellation 
term whose contribution should be added to the mass 
shift from coupled-channels to give the renormalized 
mass shift. The renormalized mass shift contains 
both a centroid correction and a fine splitting one. 
The centroid contribution will modify the quark cen- 
tral potential 0| . It is the fine splitting correction we 
will consider in this work. 
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3 Numerical calculation of fine split- 
ting 

In our calculation, the quark model is taken from 
Ref. Q . The potential parameters are: 

a s = 0.55, a = 0.175GeV 2 , m c =1.7GeV, 

C=-0.271GeV, to, = 0.33GeV, m s =0.5GeV. 

(26) 

The SHO parameter f3 is determined from the mean 
square radius of the meson state. The /3 values of 
open-charm states are 



In our calculation we take the gluon mass m g 
640MeV. This gives 



r(^(3770) ->■ DD) = 28.2MeV, 



(28) 



(3 D = 0.385GeV, (3 Ds = 0.448GeV, 



(27) 



to fit the expermental value 27. 3± l.OMeV 

To calculate the mass shift, we need further to 
know the physical mass M in Eq. (|22|) . For the char- 
monium IS", IP and 25 multiplets, we can directly 
use the experimental masses from PDG Q. For the 
2P and ID multiplets, the physical masses are the 
predicted values calculated from the assignments of 
^(3770) to l 3 D 1 and X(3872) to 2 3 P 1 . 

The mass shifts are listed in Table [2] In our cal- 
culation we take the cutoff paramter A = 800MeV. 
We also show the mass shifts without the integration 
cutoff. The cutoff reduces the mass shift by ~ 15%, 
which means that the contribution from high transfer 
momentum will be about 85% if we do not make the 
cutoff in this non-relativistic calculation. 



Table 2. The mass shifts of charmonium states in MeV. The last column lists the total mass shifts without 
the integration cutoff. 



and the f3 values of charmonium states are listed in 
Table [U 

Table 1. The f3 values of charmonium states. 



nL 


15 


25 


IP 


2P 


ID 


/3(GeV) 


0.676 


0.485 


0.514 


0.435 


0.461 



n 2S+1 Lj 


DD 


DD* 


D*D* 


D S D S 


d.d; 


did: 


total 


no cutoff 


l 3 5x 


-9 


-36 


-64 


-6 


-26 


-49 


-190 


-1359 







-52 


-47 





-39 


-36 


-175 


-1274 


1 3 P 2 


-12 


-32 


-75 


-5 


-15 


-37 


-175 


-1035 


l :i P 1 





-53 


-52 





-21 


-26 


-152 


-1021 


l 3 P 


-23 





-67 


-7 





-34 


-131 


-968 







-61 


-50 





-27 


-24 


-162 


-1021 


2 3 S 1 


-6 


-18 


-31 


-1 


-4 


-8 


-68 


-872 


2 1 S„ 





-28 


-21 





-7 


-6 


-62 


-839 


2 3 P 2 


-1 


-9 


-16 


-1 


-3 


-7 


-37 


-691 


2 :i P 1 





-17 


-10 





-4 


-4 


-35 


-716 


2 3 P 


-5 





-13 


-1 





-5 


-25 


-680 


2 1 P 1 





-18 


-10 





-5 


-4 


-36 


-701 


l 3 D 3 


-8 


-18 


-49 


-2 


-5 


-15 


-98 


-652 


l z D 2 





-40 


-33 





-9 


-11 


-93 


-665 


1 3 D 1 


-28 


-14 


-38 


-2 


-3 


-13 


-98 


-669 







-44 


-31 





-11 


-9 


-95 


-657 



The fine splittings are listed in Table [3] For 15, 
IP, 25 states, the physical mass is the experimen- 
tal mass. Then the fine splitting is calculated for 
each multiplet and listed as "splitting required" . The 
fine splitting from the quark model is calculated from 
the bare masses of the quark model which is also 



taken from Ref. p}- The fine splitting from coupled- 
channels are listed in the last column. So the to- 
tal model fine splitting is the sum of the contribu- 
tions from the quark model and from the coupled- 
channels. The results show that the calculated split- 
tings fit the "splitting required" well in 15 and 25 
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multiplets. However in the IP multiplet, the model 
splittings seem too large. 



Table 3. The physical masses and fine splittings. 



n 2S+1 Lj 


mass 


splitting 
required 


splitting 
Cj. m. 


splitting 
c. c. 


1 3 C 


QA07 


i on 


+oZ 


— 4 




/you 


— ol 


— y i 


i in 
+ LZ 


1 3 p 


Q c: c 


+ol 


+00 


1 Q 


1 3 p 


OOll 


— 10 


— iy 


+ 11 


1 3 p 


O410 


1 1 n 

— 1 1U 


i n£ 


+ol 


1 1 P 


9 £0 K 

oolo 


i n 
+0 


r 

—o 


o 

+0 


2' i S 1 


3686 


+ 12 


+14 


-2 


2 1 S„ 


3637 


-37 


-41 


+5 


2 3 P 2 


3918 


+30 


+32 


-2 


2 3 J P 1 


3872 


-17 


-17 


+0 


2 3 P 


3808 


-80 


-90 


+ 10 


2 1 P 1 


3881 


-7 


-6 


-1 


1 3 £» 3 


3798 


+6 


+8 


-2 


1 3 ^2 


3795 


+3 


-0 


+3 


1 3 A 


3773 


-19 


-17 


-2 




3793 


+0 


-0 


+ 1 



Next, we turn to the 2P and II? multiplets. This 
time, the "required spitting" is the sum of the split- 
ting from the quark model and from the coupled- 
channels. For the ID multiplet, the ^(3770) is as- 



signed to the l 3 !?! state. Then the masses of other 
states in the multiplet are calculated from the fine 
splittings as the prediction. The predicted mass of 
V-D 2 is 3793MeV. So the cc l x D 2 state is unlikely to 
be the experimental A" (3872) state even when we have 
considered the fine splitting from coupled-channels. 
So we assign the X(3872) to the 2 3 Px state and calcu- 
late the masses of the rest states in the 2P multiplet. 

4 Summary 

We have calculated the fine splitting in charmo- 
nium spectrum in the quark model with the chan- 
nel coupling effect. The open charmed meson- meson 
channels below 4GeV, including DD, DD*, D*D* 
and D s D a , D S D*, D* S D* S , are considered. The 
current-current nonlocal interacting action is con- 
structed from the color confinement and the one- 
gluon exchange interaction in the quark model. Using 
the massive gluon propagator from the lattice calcu- 
lation in the nonperturbative region, the coupling in- 
teraction is further simplified approximately to the 
four-fermion interaction. The numerical calculation 
still prefers the assignment 1 ++ of A" (3872) after we 
consider the fine splitting effect from the coupled- 
channels. The 2P and ID charmonium spectrums are 
estimated from the assignments of 1 3 D 1 to ^(3770) 
and 2 2 P 1 to X(3872). 

We would like to thank professor Shi-Lin Zhu for 
useful discussions. 
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Appendix 



For the confinement, F^ nf can be represented as 



The Partial- Wave Amplitudes 

The partial- wave amplitude is the sum of contribution 
from the confinement and from the coulomb interaction: 

a ^conf 9 *coul- 

mi mi 



pis 



In the following, 



k (/3i+/?I) fc/2: 



m q + m c 

m c 
m q + m c ' 



(29) 



(30a) 
(30b) 
(30c) 



F l c L f =—Fi{p)C h 



where C a is a spin-orbit recoupling coefficient 



(31) 



C is =(-1) 



sc + s+Ia+JA 



SA S 1 

I Ia 3a 



1 


1 


SB 


2 


2 


1 


1 




2 


2 


sc 


SA 


1 


s 



x \/6(2s + 1) {21 A + l)(2s A + l)(2s s + 1) (2s c + !)■ 

(32) 

The (p) is the polynomial of transfer momentum p: 



F P (2S^IS+IS)-- 
F S (1P-HS + 1S)-. 
F d (lP^lS+lS)- 
F S (2P^IS+IS)- 
F d (2P^lS+lS)- 

F P (1D^IS+IS)- 
F f (lD-HS+lS): 



3V3 

[ ^{ [U7D?'~SD?) + 6UD¥ - Df)] p- 2g&D? + 2^Df)p 3 } 



9 



X 2 ^D?-UicD? + 2t q Df )p>} 



9^ 
16 

9^/3 
8 



9^ 
8^/2 



9V15 
Wy/2 
45 
16 



15\/3 



(33) 
(34) 
(35) 

MtcD? + 2£ q D?) P 2 (36) 
{ 15(D 3 7 B - D\ 7 ) - 5£ c [£ c (3i?| s — -Dg 7 ) + 2^(S? 3 - Z^ 5 )] p 2 + 2£(&£>?? + 2^ q Df )p 4 } (37) 
{£c [? C (9I>9 35 ~ 5D}/) + WUD'f - DY')] P 2 - 2e c (^Dll+2^ q Dl 3 )p 4 } (38) 

^cDfp - £ {£ c Df + 2t q D? )p 3 ] (39) 

(40) 



£{^Df + 2^ q D 2 r 3 )p 3 



For the one-gluon exchange, F c * ul is further decom- 
posed to 

Fiom = —Fu(p)C ls + — F 2l (p)C ls -—F u (p)C l 2 s , (41) 

TTlq Tfbc TTlc 

where C l 2 is another spin-orbit recoupling coefficient. 
• sa = sb = sc = 1 



• SA = S B = 1, Sc =0 



c ^i_ ( 1) u +J , 4+ i v/3(2U + l) Jl 1 1 



i ^A JA 



grf. = ^Ia+ja v /2(2s + 1)(2/a + 1) 

• SA — 1, SB — SC — 



s Ja 1 



SA=0 



C'' s=1 = 0. 



2(21a + 1) 
2JA + 1 ' 



The polynomials Fu(p) and F 2 i(p) are: 
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F lp (lS -> 15 + 15) = -JL=&(£>™ - Pf )p (42) 

F 2p (15^15 + 15)=^=£ c (Pf + Pf)p (43) 

F lp (25 -»■ 15 + 15) =- ^ [£ C (7P 25 - 3P 43 + 3Pf - 3Pf> + 2£ 3 (P 43 - P?V] (44) 

F 2p (25 -> 15 + 15) = ^ [Cc(7I>? 5 - 3Pf - 3P 23 + 3Pf)p + 2£ 3 (Df + P? V] (45) 

F ls (lP^ 15 + 15)=^ [3D ls +ec{D 33_ D is )p 2 ] (46) 

F 2s (IP ^ 15 + 15) = - ^ [3^ 5 + £ 2 (Pf + Df) P 2 } (47) 

F 1(J (1P ^ 15 + 15) = - ^=£ c 2 {Df - Df)p 2 (48) 

F 2 4lP^15 + 15)=^=£ c 2 (P? 3 + P 13 )p 2 (49) 

F u (2P -S- 15 + 15) = ^= [l5(Pf - P^ 7 ) - 5£ 2 (3Pf + Pf - Pf - P^ 5 )p 2 - 2^(Pn - Df )p 4 ] (50) 

9V15 

P 2s (2P ->• 15 + 15) = -?= [l5(Pf - P 17 ) - 5£?(3Pf - Df - P? 3 + P 15 )p 2 - 2£(£>g + Df)p 4 ] (51) 
9v 15 

F ld (2P ^ 15 + 15) = - [£ 2 (9Pf - 5P^ 5 + 5Pf - 5P 9 53 )p 2 + 2£ 4 (P3 i3 - P 3 , V] (52) 
9v 15 

F 2d (2P 15 + 15) = 01 (9Pf + 5P 15 - 5Pf - 5Pf )p 2 + 2£ 4 (P?? + P 9 33 )p 4 ] (53) 

F lp (lP -> 15 + 15) =^ [5CePf P + d(^9 3 ~ £>tV] (54) 

F 2p (lP^15 + 15)=-±^ [5£ c P 2r > + £cW + .D?V] (55) 

F 1/ (lP^15 + 15)=-^=e c 3 (P 43 -P 2 V (56) 

F 2/ (1P^15 + 15)=-^ 3 (P 9 43 + P 7 2 V (57) 



